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Abstract 

In this note, we discuss a class of time-dependent Hamilton-Jacobi equations depending on a function of time, 
this function being chosen in order to keep the maximum of the solution to the constant value 0. The main 
result of the note is that the full problem has a unique classical solution. The motivation is a selection-mutation 
model which, in the limit of small diffusion, exhibits concentration on the zero level set of the solution of the 
Hamilton-Jacobi equation. The uniqueness result that we prove implies strong convergence and error estimates 
for the selection-mutation model. To cite this article: S. Mirrahimi, J.-M. Roquejoffre, C. R. Acad. Sci. Paris, 
Ser. I 340 (2015). 

Resume 

Unicite pour une classe d’equations de Hamilton-Jacobi avec contraintes. Dans cette note, on discute 
une classe d’equations de Hamilton-Jacobi dependant du temps, et d’une fonction inconnue du temps choisie pour 
que le maximum de la solution de I’equation de Hamilton-Jacobi prenne tout le temps la valeur 0. Le resultat 
principal de cette note est que le probleme complet admet une unique solution classique. La motivation est un 
modele de selection-mutation qui, dans la limite d’une diffusivite nulle, prsente une concentration sur la ligne 
de niveau 0 de la solution de I’equation de Hamilton-Jacobi. Le rsultat d’unicite que nous demontrons implique 
une convergence forte avec estimations d’erreur pour le modele de selection-mutation. Pour citer cet article : S. 
Mirrahimi, J.-M. Roquejoffre, C. R. Acad. Sci. Paris, Ser. I SjO (2015). 
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Version frangaise abregee 


On presente dans cette note un resultat d’unicite pour le probleme de Hamilton-Jacobi suivant, d’in- 
connues {u{t, x), I{t)) : 


dtu = \Vu\^ + R{x,I), {t>0,xGR‘^), maxu{t,x) = 0, 

X 

/(O) = /o > 0, ^(0, x) = ^ 0 ( 0 ;). 


( 1 ) 


La donnee R(x, I) verifie des hypothses de stricte concavite par rapport a a; et de monotonie par rapport 
a I explicites plus bas. La donnee initiale uq verifie aussi des hypotheses speciales de concavite. Ainsi, la 
fonction I{t) doit etre choisie pour que la solution u{t,x) de I’equation de Hamilton-Jacobi ait, a chaque 
instant, un maximum egal a 0. 

Nous avons alors le 

Theorem 0.1 On choisit R € et on suppose I’existence de Im > 0 tel que max R{x, Im) = 0 = 

R{0,Im)- De plus, R est supposee strictement coneave et, pour |x| grand, eomprise entre deux paraboles. 
La donne initiale uq est egalement a dcroissanee quadratique et strictement coneave. 

Le probleme (1) a une unique solution (u,I), oil u est une solution classique de I’equation de Hamilton- 
Jaeobi. De plus u e (i?+; O {R+; L'^^iR‘‘)) x lTb“(R). 

L’existence pour (1) a ete demontree en plusieurs endroits, voir par exemple [8] ou [1]. L’unicite est done 
notre resultat principal, e’etait un probleme ouvert. L’unicite etait en effet connue seulement pour un cas 
trfe particulier (voir [8]). 

Le modele (1) intervient dans la limite £ ^ 0 des solutions de 

dtUg — eAue =—R[x,Ie{t)) {t > 0, X Ie{t) = / ^p{x)ns{t,x)dx, (2) 

£ jR'i 


ou ne{t,x) est la densite d’une population caracterisee par un trait biologique x d-dimensionnel. La 
competition pour une ressource unique est representee par Rlt), ip > 0 reguliere donnee. Le terme R{x, I) 
est le taux de reproduction. Les hypothses de concavite sont des hypotheses techniques, mais pertinentes 
au plan biologique. Par une transformation de Hopf-Cole = exp (ug/e) on se ramene a I’equation sur 
Ug suivante : 

dtUg = eAug+ \Vug\'^ + R{x,Ig) (3) 


qui, dans la limite e —> 0, donne I’equation sur u. On s’attend alors a ce que Ug se concentre aux points 
ou u est proche de 0. Et, dans cette limite, Ig apparait comme une sorte de multiplicateur de Lagrange. 

La convergence de (3) vers (1) a une sous-suite pres est connue depuis [8]. Le Theoreme 0.1 donne la 
convergence de toute la famille, ainsi que des estimations d’erreur. Soit Xg(t) le point ou Ug(t, .) atteint 

son maximum. On suppose I’existence de tel que 0 < < 4(0) := [ Pj{x)n°{x)dx < Im, et on 

suppose 

with Mo G C'^(M‘’*) and maxuo(a;) = 0. 


0 _ ^ „uo/e 

'e ^ g.d/2 ’ 


Le resultat est alors le 

Theorem 0.2 Soit Ug la solution de (2) et Ug definie par (3). Nous avons les developpements asympto- 
tiques suivants 


Ig = 1 1-eh I-o{l), Xg = X I-exi + o{l), Ug 


u -I- £log(—) -I- £Mi -I- 0 ( 1 ). 
£2 


Les termes h, xi et ui vont etre presentes dans [7]. Ce resultat implique le corollaire 
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Corollary 0.3 Nous avons Vapproximation suivante pour : 

ne{t, x) = ^ (^exp(Mi + + o( 

En particular, lorsque £ —>■ 0, toute la suite (ne)e converge : 

ng{t,x) — p{t) S[x — x{t)) au sens des mesures, 


avec p{t) 


m 

4 ){x{t)) ■ 


En d’autres termes, la population se concentre sur un trait dominant qui evolue avec le temps. On note 
que la convergence de a une sous-suite pres etait deja etablie dans [5]. 

Tous ces resultats seront detailles dans [6] et [7]. 


1. Introduction 


The purpose of this note is to discuss uniqueness in the following problem, with unknowns (/ (t),u{t,x)): 


dtu = \Vu\'^ + R{x,I) (t > 0, a; G ma,xu{t,x) = 0, 

X 

/(0)=/o>0, u(0,a:) = Mo(a;), 


( 4 ) 


where Iq > 0 and uq is a concave, quadratic function: 

-Lj^-L^\x\^ <uo{x) <Lo-Li\x\^, - 2 L^ < D^uq <- 2 Li, D^uoGL^iW^). 

The constraint on the maximum of u{t, .) makes the problem nonstandard. Our main result is 
Theorem 1.1 Choose R G C^, and suppose that there is Im > 0 such that max R(x, Im) = 0 = i?(0, Im)- 

Also assume the following concavity and regularity properties for R: 


-Ki\x\^ < R{x,I) < Ko-Ki\xf, for 0 < / < Im, 

-2K^ < D'^R{x,I) < -2Ki < 0 and D^R{-,I) G for 0 < I < Im, 

-K^ < diR < -K2, 

7)1 + \d^RixiXi{x,I)\ < K 3 , forO<I< Im, and i, j = 1,2, ■ ■ ■ ,d. 

Problem (4) has a unique solution (u,I), where u solves the Hamilton-Jacobi equation in the classical 
sense. Moreover u € L'^^{R+;Wi 4 f°lw^)) O (7?+; (R'^)) x lTb°°(IR). 

Existence to (4) has been proved in various contexts (see [8,1,5]). Thus, our contribution is uniqueness, 
which has up to now been an open problem. The uniqueness has indeed been known only for a very 
particular case (see [8]). 

The rest of the note is organized as follows. In Section 2, we explain the motivation and, in particular, 
the meaning of the various assumptions. In Section 3, we revisit existence for (4), which will entail an 
unconventional ODE formulation for uniqueness. Section 4, which is the main part of the note, provides 
a fairly complete sketch of the uniqueness proof. In Section 5, we give an application. 


2. Background and motivation 

Model (4) arises in the limit £ ^ 0 of the solutions to the problem 
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( 5 ) 


dtUg — eAug =—R[x,Ie{t)) {t > 0, X Ie{t) = / ^p{x)ne{t,x)dx, 

£ jR'i 

where ng{t,x) is the density of a population characterized by a d-dimensional biological trait x. The 
population competes for a single resource, this is represented by Igit), where -ip is a given positive smooth 
function. The term R{x, I) is the reproduction rate; it is, as can be expected, very negative for large x 
and decreases as the competition increases. Such models can be derived from individual based stochastic 
processes in the limit of large populations (see [2]). The concavity assumption on i? is a technical one, 
although biologically relevant. The Hopf-Cole transformation Ug = exp (ug/e) yields the equation 

dtUg = eAug+ \Vugf + R{x,Ig) ( 6 ) 

which, in the limit £ ^ 0, yields the equation for u. Now, Ig being uniformly positive and bounded in s, 
the Hopf-Cole transformation leads to the constraint on u. Moreover, one expects that Ug concentrates at 
the points where u is close to 0 and the function Ig appears, in the limit, as a sort of Lagrange multiplier. 

This approach, based on the Hopf-Cole transformation, to study (5) has been introduced in [4] and 
then developed in different contexts (see for instance [8,1,3,5]). Long time asymptotics of such models 
have also been studied in [9] and the references therein. 


3. Existence 

Existence to a solution to (4) is obtained by letting £ —> 0 in (6). The main step is the 
Theorem 3.1 (uniform estimates for Ug, [5]) There exists Im > 0 such that 0 < Im A Ie{t) A Im + Ce^. 
Moreover we have the following estimates on Ug 

j —Lq — — e2dLit < Ug{t,x) < Lq — Li\x\^ + (i^o + 2dsLi) t, 

\L^-2tK^<D\g{t,x)<-2Li, \\D^Ug{t,-)\\L^ <C{T), for t G [0,T]. 

The bounds for Ug can be obtained for any uniformly bounded function Ig, not only for that of (5). This 
remark will be an important ingredient of the uniqueness proof. 


4. Uniqueness 

For a given continuous function 1(f) such that 0 < I{t) < Im, one may construct a solution of dtu = 
|Vw|^ -I- R{x,I) with initial datum mq. Just as in Theorem 3.1, this solution satisfies estimates (7). And 
so, u(t ,.) being strictly concave and quadratically decreasing, there exists a unique function x{t) such 
that u(t,x(t)) = maxu(t,x). Assume that Kf) is chosen such that u(t,x(t)) = 0. Then, from the equation 

on u we deduce that R{x{t),I{t)) = 0. Notice also that, because djR < 0, we have R{x(^t),0) > 0. 
Finally, differentiating 'Vu{t,x(t)) = 0 and plugging in the equation for u we obtain an ODE for x: 
x{t) = {-D‘^u{t,x{t))) ^ V,^R{x{t),I{t)). 

The idea is thus to change the constrained problem (4) by the following slightly nonstandard differential 
system: 

ri?(T(t),/(t)) = 0, forte[0,r], 

<.x{t)={-D‘^u{t,x{t))) ^ VxR(x{t),I{t)), forte [0,r], (8) 

[dtu=\yu\‘^ + Rix,I), in [0,T] X 
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with initial conditions 


7(0) = Jo, u{0,-) = uo{-), x{0) = xq, such that 7?(a;o, Jo) = 0. (9) 

Note that (8) is really a differential system because the assumptions on R imply that I{t) can implicitely 
be expressed in terms of x{t). And it is slightly nonstandard because x solves an ODE whose nonlinearity 
depends on u. Finally, note that, as soon as u satisfies the concavity and regularity estimates (7), system 
(8) is equivalent to the constrained problem (4). 

This suggests to use a simple fixed point argument to prove uniqueness to (8) (and so, to (4)). Which 
in turn suggests to set up the following scheme: starting from x{t) G (^([O,T],such that a;(0) = Xq, 
where R{xo, 0) > 0. Let J(f) solve R{x{t),I{t)) = 0 on [0, T] with i?(xo, Jq) = 0. Let v{t, .) be the unique 
solution to 

= iVup + i?(x, J), v{0,x) = uo{x). (10) 

Let y{t) solve y{t) = (—J?^v(t, x(t))) ^ V 2 ;i?(x(t), J(t)) on [0,r] with initial datum j/(0) = Xq- Setting 
y := 4>(x), we notice that uniqueness is proved as soon as we have proved that has a unique fixed point. 
One additional feature about a solution {I,u,x) of (8): 

Lemma 4.1 The function I{t) is increasing. 

We claim that our problem reduces to proving the 

Theorem 4.2 There exists C > 0 universal and 6 > 0, which is small as R(xo,0) tends to 0, such that 
^ is a contraction from (^([O, 5], Bc{xo)) to itself; here Br{a) denotes the ball of centre a G and radius 
r > 0. 

Note indeed that, by Lemma 4.1, we have, because djR < 0: 

R{x{S), 0) = R{x{S), 0) — R{x{S), I{S)) > cl{5) > cIq, 

for some universal c > 0. Hence Theorem 4.2 can be iterated to yield global existence and uniqueness. 

Let us give an overview of the proof of Theorem 4.2. For J G C ([0,5]; [0, Jm]), let V(!) be the (unique) 
solution of (10). The main step is the following 
Lemma 4.3 Jet Ji, J 2 G C ([0,5]; [0, Jm])- Then 

||F(Ji) - F(J2)||w2.~([0.5]xR‘i) < C'll-^l - -^2 ||l°°([0.<5])^- 

This lemma, once proved, opens the way to Theorem 4.2. Indeed the equation R{x, J) = 0 yields a 
smooth mapping x ^ I, and J F is a Lipschitz mapping thanks to Lemma 4.3. Moreover, the equation 
y{f) = [—D‘^v(t,x{t))) VxR{x{t),I{t)) yields a Lipschitz mapping v i-G j/ by the estimates for v given 
by Theorem 3.1. 

Lemma 4.3 is more involved. If Ji and I 2 are as in the assumptions of the lemma, the function r = 
F(Ji) — F(J 2 ) solves 

(dtr=(Vvi+Vv 2 )-Vr + R(x,l 2 )-Bix,Ii), in [0, 5] x 
'|r(0,a;) = 0, for all x G 

with Vi =V{Ii). Note that the above equation has a unique classical solution which can be computed by 
the method of characteristics. The characteristics solve 

7(t) =-Vui(t,7) - Vw2(L7), (12) 

and, due to the estimates of Theorem 3.1, they exist globally. So, one may successively express r given 
by integration along characteristics, and estimate its derivatives recursively. 
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5. Application 


Convergence of (6) to (4) had already been proved in [8,1], along subsequences. The uniqueness part of 
Theorem 1.1 yields the convergence of the full family of solutions of (6), instead of convergence along 
a subsequence. Moreover it allows an expansion of Ig, Ug and Xg (the maximum point of Ug at each time) 
in terms of e. Here are the results. 

Assume that there is such that 0 < 7° < /e(0) := / 'ip{x)ng{x)dx < Im, and that 


’’’ u°le 
C gd/2 ^ ’ 


,0 ^ ^ 


with G and maxu^(a;) = 0. 

a;GR‘^ 


The result is the 

Theorem 5.1 Let rig he the solution of (5) and Ug he defined hy (6). We have the following asymptotic 
expansions 

T 

Ig = I + ell + o{l), Xg =x-\-exi + o{l), Ue = u + elog(—) + £Mi + o(l). 

£2 

The terms /i, xi and mi will be provided in [7]. This yields the corollary 
Corollary 5.2 We have the following approximation for Ugi 

fit, x) 

exp yun^(,,xy -f - 

£2 

In particular, os £ —)• 0, the whole sequence {ng)g converges: 

ng{t,x) —>■ p{t) 6[x — x{t)'), weakly in the sense of measures, 


, , r / ! , , u(t,x)^ 

= — (exp {ui{t,x) H- -) + ( 

P2 V 6 


with pit) = , A' . 

fi(.x{t)) 

In other words, the population density concentrates on a dominant trait which evolves in time. We note 
that the convergence of Ug along subsequences was already established in [5]. 

The above results will be detailed in [6] and [7]. 
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